Extracting Boer-Mulders functions from p + D Drell-Yan processes 
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We extract the Boer-Mulders functions of valence and sea quarks in the proton from unpolarized p + D Drell- 
Yan data measured by the FNAL E866 Collaboration. Using these Boer-Mulders functions, we calculate the 
cos 20 asymmetries in unpolarized pp Drell-Yan processes, both for the FNAL E866/NuSea and the BNL Rel- 
ativistic Heavy Ion Collider (RHIC) experiments. We also estimate the cos 2<p asymmetries in the unpolarized 
pp Drell-Yan processes at GSI. 
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I. INTRODUCTION 



Recently the E866/NuSea Collaboration at FNAL has mea- 
sured llJ] the cos 20 angular distribution of Drell-Yan dimuons 
in p + d interaction at 800 Ge V/c. The magnitude of the asym- 
metry turned out to be about several percent. This cos 20 an- 
gular dependence, together with other angular correlations in 
Drell-Yan processes, constitute remaining challenges which 
need to be understood from QCD dynamics. Indeed, even 
before these nucleon-nucleon interaction measurements, the 
NA10 |2] and E165 |3] Collaborations had also measured the 
cos 20 angular dependence in n + N Drell-Yan processes. The 
magnitude of the angular dependence is around 30% for mod- 
est transverse momenta of the lepton pair, a result which vi- 
olates the so-called Lam-Tung relation [4], predicted by per- 
turbative QCD. Several attempts have been made to interpret 
these data, including QCD vacuum effects (J, |6(] and higher- 
twist mechanisms (71 |8fl. Furthermore, in the last decade a 
significant efforts have been put forward on the understand- 
ing of the cos 2(f> angular dependence from the view point 
of the transverse momentum dependent (TMD) Boer-Mulders 
function hf [90. In fact, it was shown [10] by Boer that 
the angular dependence can be related to the product of two 
functions hf, each of which comes from one of the inci- 
dent hadrons. These hf functions describe a correlation be- 
tween the transverse spin and the transverse momentum of 
a quark inside an unpolarized hadron, which originates from 
inital/final state interactions lflll[l"2tl . and which in turn is re- 
lated to t he gaug e invariance of the TMD distribution func- 
LL4j,[15I]. Another distribution having the same QCD 



tions _ . 

origin of hf is the Sivers function JT6fl, which plays a impor- 
tant role in the single spin asymmetries (SSA) 11711 measured 
in semi-inclusive deep inelastic scattering. Measurements of 
Sivers SSA at HERMES H B and COMPASS H 
have been used in order to extract the Sivers functions by sev- 
eral groups I22L I23L l24l I25L 12611 . All these fittings show that 
the Sivers functions for u and d quarks have similar size but 
opposite sign. 



In contrast to the Sivers function case, detailed knowledge 
of hf is more difficult. The reason is that in processes in 
which the Boer-Mulders function contributes, one finds that 
it is always convoluted with itself or with other chiral-odd 
functions, such as the transversity or the Collins fragmenta- 
tion function 12711 . Despite this, some theoretical calculations 
and phenomenological analysis d M, M, M, M, HI 
l34l l35l [36l H2I [Hfl have been performed on hf. An inter- 
esting problem is its flavor dependence [39]. Recently lattice 
calculations B4011 and ap proa ches based on generalized par- 
ton distributions (GPD) HH |H] suggest that the sign of hf 
for both u and d quarks and the corresponding magnitudes 
are of similar sizes (See also the calculation in Ref. 13710 . 
In this paper, we will extract the Boer-Mulders functions 
from the unpolarized p + d Drell-Yan data measured by the 
E866 collaboration, with the assumption that the cos 2<p an- 
gular dependence is produced solely by hf. We use the sim- 
ple parameterization for the Boer-Mulders function given by 
hf- q (x, p]) = H q x c (\-x) /*(*) exp (-pI/pIJ to fit the data 
points of the asymmetry coefficient v vs pj, x\ and X2, where 
q stands for the u, d, u, d quarks. The data points of v versus 
m mi - Q 2 an d x f provide then a cross check of our fitting. Us- 
ing the resulting valence and sea quarks Boer-Mulders func- 
tions, we estimate the cos 20 asymmetries in unpolarized pp 
Drell-Yan processes for both FNAL E866/NuSea and RHIC, 
and we also give predictions for the cos 20 asymmetries in 
unpolarized pp Drell-Yan processes at GSI. 



II. EXTRACTING THE BOER-MULDERS FUNCTION 
FROM UNPOLARIZED p + D DRELL-YAN DATA 

The angular differential cross section for unpolarized Drell- 
Yan processes has the general form 
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where 8 and <f> are, respectively, the polar angle and the az- 
imuthal angle of dileptons in the Collins-Soper frame l43ll . 
The coefficients A, p. and v do not depend on these angles, and 
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TABLE I: Best fit values of the Boer-Mulders functions 



H u 


3.99 


H d 


3.83 


H-u 


0.91 


Hi 


-0.96 




0.161 


C 


0.45 


rld.o.f. 


0.79 



for scattering that has azimuthal symmetry their values are 
H = v = 0. 

This angular distribution has been measured in muon pair 
production by pion-nucleon collisions: n~N — > /i 4 pTX, with 
N denoting a nucleon in deuterium or tungsten, and for a 
iC beam with energies of 140, 194, 286 GeV Jl and 252 
GeV H. The experimental data show large values of v, near 
30%. The most recent measurements of the angular distribu- 
tion were performed by the E866 Collaboration JTt], in p + d 
Drell-Yan processes at 800 GeV/c. The measured v is about 
several percent, a result which can not be explained by pertur- 
bative QCD. As proposed by Boer [ 10], the non-zero cos2</> 
term can be produced by the product of two hfs, each coming 
from one of the two incident hadrons. In the case of the p + D 
Drell-Yan processes, the coefficient v can be expressed as lf3ltl 



2T[ X (e\h\' u + e 2 X' d )(hf + hf)] + (q «-> q) 

VpD = = 

nielft + e 2 ff)(ff + ff )] + {q ~ q) 
where we have used the notation 

T[- ■■] = f d 2 p±d 2 k ± 6 2 ( Pl _ + k ± - q ± ) x {• ■ ■}, 

and 

X(P±, k±) = (2h P± h k ± - p ± kJ/M 2 , 



(2) 



(3) 



(4) 



where h = qr/Qr and M is the mass of the nucleon. The in- 
teracting quarks coming from the incident hadrons have trans- 
verse momenta p ± and k ± , and Qj is the absolute value of the 
transverse momentum of the photon q^. To arrive at Eq.© 
we have used the isospin relation: 



(5) 



for f\ and hf. From this point of view, the angular distribu- 
tion coefficient v measured by the E866 Collaboration pro- 
vides an opportunity to extract the Boer-Mulders function of 
the nucleon. For this purpose we parameterize hf(x, p^.) in a 
factorized form as 



^• ? (*,Pr)=^i M 



exp(-p 2 T /p 2 h J 
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(6) 



which is a Gaussian model for the transverse momentum de- 
pendence of the Boer-Mulders functions, with width p 2 bm . As 



in Ref. [10], the x-dependence of the Boer-Mulders functions 
are modeled relating its behavior with that of the unpolarized 
distribution functions, as 



h\<\x) 


= H u x c (l 
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h\ 4 {x) 


= H 3 x c {\ 


-x)ff(x), 


(10) 



where f y q {x) is the well-known unpolarized integrated distri- 
bution function. Therefore we have H u , Hd, H[i,H d ,p 2 }m and c 
as parameters in our parametrization. The coefficient (1 —x) is 
included in order to give the correct large-x behavior J44ll for 
hf compared with the unpolarized distribution. The small-x 
behavior of hf compared with that of f\(x) is modeled as x c . 

The TMD unpolarized distribution function f\ (x, p^) is also 
given in a Gaussian form 



y?u.i»r) = y?w 



exp (-p 2 T /p 2 unp ) 



(11) 



From Eqs. © to (HU, we can de-convolute the transverse 
momentum integration in Eq. (|2| and arrive at 



V pD (x U X2, Qt) 



p 2 unp [a]Q 2 T exp(-^ 
2MVL^]ex P (- 5 | ; )' 



(12) 



where 



[a] = x\(\ - Xl )x c 2 (l - X2)(4H u f?( Xl ) + H d ff( Xl )) 
(H a f, u (x 2 ) + H s f«(x 2 )) + (q~ q), 

m = (4/i"(xi) + ff (xi))(/ffe) + ff(x 2 )) 
+(q <-» q) 



(13) 



(14) 



Therefore the coefficient v vs Qj can be obtained from (fT2l) 
by integrating the numerator and the denominator over x\ and 
x 2 over, respectively: 

PlnpJ dx \ J dx 2 [a]Q 2 T expl-^r) 
v p d(Qt) = T^TT-v (15) 

2M 2 P 2 m fdx 1 fdx 2 m^P (-|y 

In this equation, x\, x 2 and Qj satisfy the relation x\x 2 s = 
Q 2 + Qj, where s and Q are the cm energy squared and the 
invariant mass of the lepton pair respectively. In a similar way 
we can obtain expressions for v p d(x\) and v p d(x 2 ). 

The E866/NuSea Collaboration has measured v p d vs Qt, 
x\, x 2 , Xf and m m in the following kinematical region: 

4.5 GeV < Q < 9 GeV and 10.7 GeV < Q < 15 GeV, 
0.1 < x, < 0.9, 0.02 <x 2 < 0.24 

We therefore use the above expressions for v p d(Qj), v p d{x\) 
and v p d(x 2 ) in order to fit the experimental data of v vs Qj, 
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X[ and X2 in pD processes measured by E866/NuSea Collab- 
oration. We perform the fitting using the MINUIT program. 
For the unpolarized distribution function f^(x) we adopt the 
MRST2001 (LO set) parametrization with QCD evolu- 
tion taken into account. We use p 2 np = 0.25, following the 
choice in Ref. l46ll . a value which was obtained by fitting the 
azimuthal dependence of the SIDIS unpolarized cross section. 
Notice that these values are assumed to be constant and flavor 
independent. 

The best fitting values for the parameters are given in Ta- 
ble U and in Fig. [TJwe show the fitting result compared with 
the E866/NuSea data. Notice that the results for the Boer- 
Mulders functions given in Table U are obtained for the Drell- 
Yan process, and to obtain the corresponding Boer-Mulders 
functions in SIDIS one should reverse their signs II 1 511 - There- 
fore our results agree with the expectation that the Boer- 
Mulders functions for u and d quarks in SIDIS are negative 
and have the same sign. We also want to emphasize that 
the Boer-Mulders functions we extracted are within a positive 
bound [471 



\p T hf(x, )| 



M 



i fi(x, Pr)- 



We also find that at 



Pt 



ylpl P P 2 bm /(2( P 2 mp -p 2 bm )) = 0.45GeV, 



(16) 



(17) 



the ratio prhf(x, pj.)/(M/i(x, p^.)) has its maximum for all 
x. As an example, in Fig. [2] we show a comparison of 
prxhf(x,Pj)/M and xf\(x, p^.), for both u and d quarks at 
Pt = 0.45 GeV and Q — 1 GeV, and therefore we conclude 
that our results obey the positive bound for all pj. 

As a cross check, we give results for the coefficient v versus 
Xf and m w in pD Drell-Yan processes for E866/NuSea, us- 
ing the extracted Boer-Mulders functions and comparing them 
with data, as shown in Fig. [3] In the calculation we use the re- 
lation xi/2 = (+Xf + (Q 2 + Qj)/s)/2. The results in the figure 
show good agreement with data. 



m. PREDICTIONS FOR pp AND pp PROCESS 

In the previous section we have performed a fitting on the 
unpolarized pd Drell-Yan data, and obtained a set of Boer- 
Mulders functions. In this section, we will use the extracted 
functions to predict the result in other processes. We first 
focus on unpolarized pp Drell-Yan processes. The cos 20 
angular distribution in this process can be also measured at 
E866/FNAL, and the expression that we get for v using our 
parametrization of the Boer-Mulders functions is 



,,(xi,x 2 ,Q T ) = 



Pimp 



M(g«p(-£) 



2MV te [A]«p(-£)' 



(18) 



+H d ft(x 1 )H i f((x 2 )) + (q~ q), 
[A] = (4f?( Xl )f?(x 2 ) + ff(x l )ff(x 2 )) 
+(q <-> q) 



(19) 



With the Boer-Mulders functions for both valence and sea 
quarks in the proton that we obtained above, and taking the 
same integration regions for Q, x\,x 2 in evaluating Eq. ( fT8l , 
we get the results for cos 20 asymmetry in unpolarized pp 
and pD Drell-Yan processes that are shown in Fig. [4] (solid 
and dashed lines respectively). 

The results in these two processes present little difference 
between them. 

Another promising testing ground of the cos 20 asymmetry 
in unpolarized pN Drell-Yan processes is BNL RHIC B48I1 . 
since this experiment is also feasible at this accelerator. With 
the Boer-Mulders functions given in Table U we also esti- 
mate the cos 20 asymmetry v at RHIC for an unpolarized ex- 
periment, with kinematical constraints t/s = 200 GeV and 
-1 < y < 2. Here y is the rapidity defined as y — lln—, with 

xi/2 - ^j Q + ^ T e±y ■ After integrating over X\ in Eq.([T8]i, we 
get the Q r -dependent asymmetries for 2 = 4 GeV (solid line) 
and Q — 20 GeV (dashed line) shown in Fig. [5] 

We should notice that in the extraction of the Boer-Mulders 
functions with the MINUIT program, the coefficients of Boer- 
Mulders functions for valence quarks H u and Hd are always 
coupled with the coefficients of sea quarks H a and H d , as seen 
in Eq. ( fT3l . We cannot separate them except by introducing 
free coefficients co q with possible flavor dependence. 

With the free coefficient co u — <dj = to, the Boer-Mulders 
functions are modeled comparing its behavior with that of the 
unpolarized distribution functions, as 



h\' u {x) = 


ojH u x c (\ ■ 




(20) 


h\ 4 {x) = 


coH d x c (\ ■ 


-x)ff(x), 


(21) 


hf{x) = 


-H- U x c (l 


-x)fl'(x), 


(22) 




to 




h\' d (x) = 


-H d x c {\ 

CO 


-x)ff(x), 


(23) 



The results we get above for the v p d,v pp are independent 
of the free coefficient co because when we use Eqs. d20l-(l23l) 
into the v p d, v pp in Eqs. ( fT2l and ( 11 St . the free coefficient to 
cancels in the product of sea and valence components. 

Recently the cos 20 azimuthal asymmetry of pp Drell-Yan 
processes has also received much attention, since there have 
been proposals to study spin phenomena in polarized and 
unpolarized pp scattering at the High-Energy Storage Ring 
(HESR) of GSI BHHi. The p T , x\ and x 2 -dependent cos 20 
asymmetry v pp , obtained using our parametrization for the 
Boer-Mulders functions with the free coefficient a> are 



pp 



(x\,x 2 , Qt) 



pL,[d6 2 r exp(-^-) 
2MVLMexp(-^)' 



(24) 



where 



where 



M = x\{\ - Xl )x c 2 (l - x 2 )(4H u f?(x 1 )H a fl'(x 2 ) 



M = co 2 ^! - Xl )x c 2 (l - x 2 ){W u f u x {x,)H u ft(x 2 ) 
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FIG. 1: Fits to the pj,x\, X2-dependent cos 2cf> asymmetries v pD for Drell-Yan processes. Data are from the FNAL E866/NuSea collaboration. 
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FIG. 4: The p T -dependent cos 20 asymmetries v in both pp 
(dotted curve)and pD(solid curve) Drell-Yan processes at FNAL 
E866/NuSea, calculated with the fitted Boer-Mulders functions pre- 
sented in Table |l] 



FIG. 2: Comparison of \p T xhf(x, p|)|/M and xf\{x, pj.) for u and d 
quarks at p T = 0.45 GeV and 2=1 GeV. Here f t is a combination 
of valence and sea quark distributions. 
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FIG. 3: The g-dependent cos 20 asymmetry v pD for Drell-Yan pro- 
cesses at FNAL E866/NuSea, presenting both the experimental data 
and the results we estimate with the best fit values of the Boer- 
Mulders functions in Table [I] (line). 



Q=4GeV 
Q=20GeV 



FIG. 5: The y-dependent cos 20 asymmetries v for pp Drell-Yan 
process at RHIC, with the kinematical conditions *>fs = 200 GeV, 
-1 < y < 2, and calculated with the Boer-Mulders functions that we 
got. The dotted and dashed lines show the asymmetries at Q = 4 and 
20 GeV, respectively. 
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FIG. 6: The -dependent cos 2<f> asymmetries v for pp Drell-Yan 
process at GSI, with the kinematics: cm. energy s = 45 GeV 2 , and 
invariance mass square of the lepton pair Q 2 = 2.5 GeV 2 . The solid, 
dashed and dotted curves correspond to the free coefficient o)=l, 0.5, 
0.25 respectively. 



+H d f?(x 1 )H d jf(x 2 )) + (q~ q), 

[T] = (4/? ( Xl )mX2) + ff{Xl)ft{x 2 )) + {q~q) (25) 

In the expression above, we use charge-conjugation rela- 



tions such as: 

julP _ ^u/p jd/P _ j4/P 

f u/P =f u/P t f x iv = f m_ (26) 

In Eq.d25l), we can see that the result for this pp process de- 
pends on the free coefficient co. Therefore we give predictions 
for v P p at GSI, as shown in Fig. [6j for three different values 
of co. The kinematics in GSI can be chosen as cm. energy 
s = 45 GeV 2 and Q 2 = 2.5 GeV 2 , where the x is in the va- 
lence region and the asymmetry can be larger than the asym- 
metry in pp processes at E866/NuSea and RHIC. The coming 
experimental data can fix co to give the exact Boer-Mulders 
functions for both quarks and antiquarks. 



IV. CONCLUSION 

In summary, we have extracted the Boer-Mulders functions 
of valence and sea quarks inside the proton from unpolarized 
p + D Drell-Yan data at 800 Gev/c, measured by the FNAL 
E866/NuSea Collaboration. With the Boer-Mulders functions 
that we get, we estimated the cos let asymmetries in unpolar- 
ized pp Drell-Yan processes at E866/NuSea and RHIC. We 
also presented an estimation of the cos 20 azimuthal asym- 
metry in pp Drell-Yan processes at GSI. We hope that future 
measurement for the cos 2cf> asymmetry at pp and pp, as well 
as in SIDIS, can help to pin down the Boer-Mulders functions 
for the valence and sea quarks. 
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